Stochastic volatility and leverage effect 
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We prove that a wide class of correlated stochastic volatility models exactly measure an empirical 
fact in which past returns are anticorrelated with future volatilities: the so-called "leverage effect" . 
This quantitative measure allows us to fully estimate all parameters involved and it will entail a 
deeper study on correlated stochastic volatility models with practical applications on option pricing 
and risk management. 
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The multiplicative diffusion process known as the ge- 
ometric Brownian motion (GBM) has been widely ac- 
cepted as one of the most universal models for specu- 
lative markets. The model, started out by Bachelier in 
1900 as an ordinary random walk and redefined in its final 
version by Osborne in 1959 [0, presupposes a constant 
"volatility" a which is equivalent to a constant diffusion 
coefficient D = a 2 . However, and especially after the 
1987 crash, there seems to be ample empirical evidence, 
given by the so-called "stylized facts" , that the assump- 
tion of constant volatility does not properly account for 
important features of markets. Some relevant examples 
of stylized facts are the "leverage" and "smile" effects 
and the skewness and "fat tails" in probability distribu- 
tions 

One of the main ideas that has come out to explain 
those features is that the market dynamics has intrin- 
sically changed in the sense that volatility is no longer 
constant. It is not a function of time either (as might be 
inferred by the evidence of nonstationarity in financial 
time series ||) but a random variable. In its most general 
form one therefore assumes that the volatility a is a func- 
tion of a random process Y(t), i.e., a(t) — <p(Y(t)). We 
could make an analogy from physics saying that specula- 
tive prices S(t) evolve in a "random medium" determined 
by a random diffusion coefficient. Most of the stochas- 
tic volatility (SV) models presented up to date suppose 
that Y(t) is itself a diffusion process that may or may 
not be correlated with price and different models differ 
from each other basically in the form of the function <p 

§§• 

As mentioned above, the hypothesis of stochastic 
volatility (SV) has been suggested to explain significant 
peculiarities observed in real markets. Thus, the smile ef- 
fect, related to the implicit volatility in option prices, has 
been thoroughly studied both qualitatively and quanti- 
tatively ||, [jj. Other features such as leverage and long 
tails are less studied and some works on SV only address 
to them from a qualitative point of view [^j while oth- 
ers measure them by giving numerical coefficients, based 
on ARCH-GARCH models, for kurtosis and skewness 
However, the time evolution and structure of the lever- 



age and tails have never been investigated. Our main 
objective here is to prove that the leverage effect can be 
quantitatively explained by a wide class of correlated SV 
models. This will allow us to overturn the main objec- 
tion against SV models: the impossibility of fitting all 
parameters appearing in these models ^] which, in 
turn, opens the door to simulations of real markets with 
far reaching practical consequences on option pricing and 
risk management. 

We recall that the leverage effect has its origin in the 
observation that volatility seems to be negatively cor- 
related with stock returns which, in continuous time fi- 
nance and in terms of speculative prices S(t), are defined 
by R(t) = ln[S(i)/S(0)]. The first explanation to this 
empirical fact was given by Black |hJ and Christie |ll[] 
in the sense that negative returns increase financial lever- 
age which extend the risk of the company and therefore 
its volatility. Hence the name of "leverage effect" . Nev- 
ertheless, the cause of this effect is still unclear since an- 
other explanation is just the contrary, that is, an increase 
of volatility makes the stock riskier which produces a fall 
of demand and the price drops . 

In a very recent paper, Bouchaud et al. |Q have per- 
formed a complete empirical study of the leverage ef- 
fect for both individual stocks and indices using daily 
data. The volatility-return (negative) correlation is 
clearly shown to have a definite direction in time - a very 
confusing fact in the literature - since correlations are 
shown to be between future volatilities and past returns. 
Bouchaud et al. conclusively prove from data that the 
negative correlation decays exponentially in time, faster 
for indices than for individual stocks p2[ . 

In this letter, we present a theoretical study on these 
correlations and show that a wide class of SV models 
completely explain the leverage effect, both qualitatively 
and quantitatively, in complete agreement with exper- 
imental observations. The starting point is the GBM 
model: 

dR = fj,dt + a(t)dWi(t), (1) 
where fi is the drift and o~(t) = o~(Y(t)) is a random 
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volatility and Y(t) is a diffusion process: 
dY = f(Y)dt + g(Y)dW 2 (t). 



(2) 



In these equations Wi(t) (£ = 1,2) are Wiener processes, 
i.e., dWi(t) = £i(t)dt where are zero-mean Gaussian 
white noises with cross-correlation given by 



(3) 



(— 1 < p < 1). As is common in finance, Eqs. (|l|)-(||) 
are interpreted in the sense of Ito and for the rest of the 
paper we will follow Ito convention Jl3| . 

Bouchaud et al. jl2| , quantify the leverage effect by 
means of a leverage correlation function defined by 
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where 



£(t) = -([dX{t + T )] 2 dX(t)) 



X(t) = R(t) - pi 



(4) 



(5) 



is the zero-mean return and Z = [(dX (t) 2 )] 2 is a conve- 
nient normalization coefficient. Bouchaud et al. have an- 
alyzed a large amount of daily relative changes for either 
market indices and stock share prices and find that (L2| 



C(r) = 



-Ae~ bT 



0. 



if r > 0; 
if t < 0; 



(6) 



(A, b > 0). Hence, there is a negative correlation with 
an exponential time decay between future volatility and 
past returns changes but no correlation is found between 
past volatility and future price changes. In this way, they 
provide a sort of causality to the leverage effect which, 
to our knowledge, has never been previously mentioned 
in the literature ||, [j| . 

Let us sketch how correlated SV models are able to 
exactly reproduce this result. We first combine Eqs. (|l|) 
and (|) to produce £(r) = (a(t)dWi (t)a(t + rfdW^t + 
t) 2 )/Z. Note that when r < Ito's rules tell us that 
dW\(t) is uncorrelated with the rest of terms then, re- 
calling that (dWi(i)) = 0, we have £(t) = if t < 0. On 
the other hand, when r > 0, dW\(t + r) is uncorrelated 
with the rest and, since (dWi(t + t)) 2 = dt, we conclude 
that 



C(r) = e(T){a{t)dW 1 {t)a{t- 



2 }dt/Z, 



(7) 



where 9{t) is the Heaviside step function and Z = 
[(a 2 (t))] 2 dt 2 . Note that we have proved the existence of 
correlations between future volatilities and past returns 
but not vice-versa. Note also that this have been proved 
independently of the underlying volatility process Y(t) 
(which needs not to be a diffusion process) and of the 
specific form of a in terms of Y . 

Suppose now that Y(t) is a diffusion process given by 
Eq. (0) . As is well known any pair of correlated Wiener 



process, such as Wi(t) and Wzit), satisfy the identity 
dWi(t) = pdW 2 (t) + y/l - p 2 dW(t), where W(t) is a 
Wiener process independent of Wz(t) (therefore, W(t) is 
independent of a) . Substituting this identity into Eq. (0) 
we get C(t) in terms of the average (a(t)a(t + T) 2 £ 2 (i)). 
This average can be calculated by means of Novikov's 
theorem [^4j with the result fl5| 

2p6(r) / , , , . i , ,SY(t + r) 
C{t) = H 9 ,) ' ( a(t)a(t + r)a'(t + r) 



(8) 

where a' = da{Y)/dY and SY(t + t) /<*&(*) is the func- 
tional derivative of Y(t + r, [£2]) with respect to £2^) 
& 
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FIG. 1: The leverage effect in the Dow- Jones daily index. We 
plot the leverage function C(t) for the Dow- Jones index from 
1900 until 2000. We see that there exists a non-negligible 
correlation when r > and negligible when r < 0. Observe 
that correlation strongly fluctuates when —3 < r < 2. We 
also plot a fit with the OU SV leverage function (|TJ]) that 
helps us to estimate a and p. 

There is a wide consensus on volatility being "mean 
reverting". This means that there exists a normal level 
of volatility jTJ| to which volatility will eventually return 
||. For a general SV model such as the pres- 

ence of mean-reversion implies restrictions on the form 
of the drift coefficient f(Y). In order to include this ex- 
perimental fact in the model, the simplest choice is to 
assume that f(Y) is a linear function. That is, 



Y 



-a(Y-m)+g{Y)b(t), 



(9) 



where a < 0. The formal solution to this equation in the 
stationary state reads 



Yit) 

from which we get |l 
SY(ti 



,-a{t-t'\ 



g(Y(t'))&(t')dt', 



sut) 



= 8(T)e- aT g(Y(t))cxp 



t + T 



g'(Y(s))&(s)ds 
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Substituting this expression into Eq. (H) yields 



where 



B(r) 



C{t) = P 9(t)B(t) 



2 (a(t)a{t + r)a'(t + r)G(t, t + r)) 



and 



G(t,t + T)=g(Y(t))e X p 



t + T 



(10) 



(11) 



(12) 



In consequence, any SV model of the form given by Eqs. 
(|l|) and (||) and whose function B(t) does not increase 
faster than e aT as r — > oo, satisfies an exponentially de- 
caying leverage as expressed by Eq. @. Moreover, if 
cr(Y) is an increasing function of Y with definite sign 
and g{Y) is positive definite (or a is decreasing and g is 
negative) we see from Eq. ([l0]) that the correlation coef- 
ficient p must be negative and driving noises W\(t) and 
Wb(t) are anticorrelated. Eqs. (|lO|)-([l^) constitute the 
main result of the paper. 

The exact form of C{t) will depend on the expres- 
sion of B{j) which in turn will depend on the SV model 
chosen. Within diffusion theory, as is the case of Eq. 
(^), there are basically three different SV models ||: 1) 
The Ornstein-Uhlenbeck (OU) model where a — Y and 
g(Y) = k (a positive constant) p7| , 2) the exponential 
Ornstein-Uhlenbeck (expOU) model where a — e Y and 
g(Y) = k U and 3) the Cox-Ingersoll-Ross (CIR) model 
where a = \/Y and g{Y) = k\/Y [Q. For all these models 
the leverage function has the form given by Eq. ©. In 
the OU model and in the CIR model (the latter with zero 
mean- reversion, i.e., to = k 2 /4a) the leverage function 
£(r) is respectively given by Ju| 



£i(t) = 2fcp6»(r) 



m 2 + (k 2 /2a)e~ 



(m 2 + fc 2 /2a) 2 



pa 



-Ca(r) = 4^0(r)e- QT , 
while for the expOU model we have 



(13) 



(14) 



£ 2 (t) = 2 j ofc6'(r)exp [- 



fc 2 (e" 



- 3/4)/a] 

(15) 

In Fig. U we show the leverage effect for the Dow- Jones 
daily index (1900-2000) and plot the leverage function 
for the OU model. 

Any market model, besides being able to reproduce 
the market dynamics, must provide a systematic way of 
evaluating its parameters. Almost all current SV mod- 
els have four parameters to estimate: p, k, m, and a. 
To our knowledge, all works on SV models presented up 
till now are able to evaluate only two of them. Thus, 
for instance, Fouque et al. M estimate k and m from 



the empirical second and fourth moment of daily data 
but they can't give a clear estimation of a and p. This 
constitutes the main criticism to SV models, say, their 
inability to estimate all parameters involved. This sit- 
uation changes completely when one measures leverage. 
Indeed, k and m are obtained, as usual, from the em- 
pirical second and fourth moment. Next, by adjusting 
e~ aT to leverage empirical data we estimate a. More- 
over, comparing the theoretical and empirical leverage 
at t = + , £(0 + ), we finally obtain p. Following this 
procedure, for the Dow-Jones daily index and using the 



OU model, we have estimated k = 
to = 18.9% year" 1 / 2 , a = 0.05 day~ 



1.4 x 1Q- A days~ 



and p 



-0.58. 
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FIG. 2: Dow- Jones historical time-series (1900-2000) and 
path simulations. We show a Dow- Jones daily returns sam- 
ple path (top), its simulation by means of the OU SV pro- 
cess (middle), and the geometric Brownian process simulation 
(bottom). All parameters of the simulations are estimated 
from the whole Dow- Jones historical time-series from 1900 to 
2000. The dynamics is traced over 300 days (approximately 
one trading year) which, for the empirical path, nearly corre- 
sponds to 1999. 

As an illustration we have simulated, using Eqs. (jl]) 
and (^), the OU resulting process with the parameters 
estimated above. We follow the random dynamics of the 
daily changes of the zero-mean return X(t) and com- 
pare it with the empirical Dow-Jones time series during 
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one trading year. We have also simulated the geometric 
Brownian motion assuming a constant volatility a whose 
value is directly estimated from Dow- Jones one-century 
data. We present these results in Fig. || and observe that 
GBM cannot describe either the largest or the smallest 
fluctuations of daily returns. We nonetheless see in the 
figure that the SV model chosen describes periods of high 
volatility together with periods of very low volatility, re- 
sulting in a more similar trajectory to the Dow- Jones in- 
dex than that of the GBM. This is quite remarkable, be- 
cause we have simulated last year trajectory using all past 
data (one century) of the Dow-Jones index thus showing 
the stability of parameters. 

These results may encourage to deepen in the study 
of statistical properties of SV models showing leverage. 
Several models have been presented in the literature 
without being able to discern which is the more realistic 
one. Now, thanks to leverage correlation, it is possible 
to estimate all parameters involved in SV models. This 
will allow us to confront in detail different models with 
the empirical statistical properties of markets. Finally, 
a better knowledge of SV models has non trivial conse- 
quences on option pricing (since classical Black-Scholes 
method is still suitable in SV models [[| 0) and, more 
generally, on risk management. 
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